Collision statistics of driven granular materials 
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We present an experimental investigation of the statistical properties of spherical granular particles 
on an inclined plane that are excited by an oscillating side- wall. The data is obtained by high-speed 
imaging and particle tracking techniques. We identify all particles in the system and link their 
positions to form trajectories over long times. Thus, we identify particle collisions to measure the 
effective coefficient of restitution and find a broad distribution of values for the same impact angles. 
We find that the energy inelasticity can take on values greater than one, which implies that the 
rotational degrees play an important role in energy transfer. We also measure the distance and 
the time between collision events in order to directly determine the distribution of path lengths 
and the free times. These distributions are shown to deviate from expected theoretical forms for 
elastic spheres, demonstrating the inherent clustering in this system. We describe the data with a 
two-parameter fitting function and use it to calculated the mean free path and collision time. We 
find that the ratio of these values is consistent with the average velocity. The velocity distribution 
are observed to be strongly non-Gaussian and do not demonstrate any apparent universal behavior. 
We report the scaling of the second moment, which corresponds to the granular temperature, and 
higher order moments as a function of distance from the driving wall. Additionally, we measure 
long time correlation functions in both space and in the velocities to probe diffusion in a dissipative 
gas. 

PACS numbers: 81.05.Rm, 05.20.Dd, 45.05.+X 45.70.Mg, 45.70.-n, 51.10.+y 
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I. INTRODUCTION 

Granular material represent a type of matter not well 
defined by conventional means. Although each granular 
particle is obviously solid, an assemblage of these parti- 
cles show distinctly non-solid behavior when subjected to 
external forces 0. In the rapid flow regime, the interac- 
tion between the grains is coUisional and the system re- 
sembles a dense granular gas. Indeed, the kinetic theory 
for dense gases formulated by Chapman and Enskog ||] 
have been modified to include the dissipative nature of 
the collisions |4j. However, a number of approxima- 
tions have to be made in any calculation that can be only 
validated by experiments. Furthermore, even if key as- 
sumptions such as equipartition breakdown |[ 0|i it 
is important to have a measure of the failure to guide 
further development. 

Energy has to be constantly supplied from an exter- 
nal source to observe a steady state in granular gas sys- 
tems. Therefore, model experiments consist of granu- 
lar particles inside a container where energy is continu- 
ously injected at a side-wall M, 0, 10 . Thus gradients 



are present in experimental granular systems, which im- 
plies that care must be taken when comparing results 
to non-equilibrium kinetic theory |jri|, |l^, ^ . With ad- 
vances in high speed image acquisition, it is now possi- 
ble to obtain positions of particles several times between 
collisions. However, particle positions and velocities can 
be obtained accurately only in two dimensions by direct 
imaging thus forcing certain constraints on the geometry 
of the system. 

One of the first such experiments to investigate veloc- 
ity distribution functions (VDFs) utilized an apparatus 
in which particles are vibrated vertically inside a narrow 



transparent box ^, [Tj, |T5j. Maxwellian statistics were 
reported for the vertical and horizontal velocity compo- 
nents of the particles parallel to the plane of the trans- 
parent side-walls. Additional interactions in these sys- 
tems arise due to collisions between particles and the 
side- walls Following this work, Wildman et al. 
were able to do long time particle tracking to measure 
diffusion constants by interpreting mean square displace- 
ment data over a very broad range of density. More re- 
cently, in a similar apparatus, Rouyer and Menon [T^ 
report that their VDFs have a universal form that can 
be parameterized by a single variable, the granular tem- 
perature. A different method of energy injection utilizes 
large flat container that is vibrated vertically to excite 
a sub mono-layer of particles |l8|, |l9) . The velocity 
of the particles in the horizontal plane are measured and 
are found to follow a non-Gaussian distribution. How- 
ever, the impact of the velocity gradient in the vertical 
direction on the observed distributions are not taken in to 
account because these components cannot be measured. 

Our experiment is a variation of the vertically vibrated 
apparatus. Spherical particles are constrained to roll 
on an inclined two dimensional surface. This geome- 
try allows for a direct investigation of the interplay be- 
tween energy injected at the side-wall and the dissipa- 
tion through inelastic collisions. In addition, the inclina- 
tion reduces the effects of gravity, therefore minimizing 
shock waves. This system has been used to demonstrate 
clustering and collapse when the inter-particle collision 
frequency is much greater than particle-driving wall col- 
lision frequency Recent works have explored a full 
range VDFs, from very near Gaussian behavior to highly 
non-Gaussian distribution functions, as well as velocity 
correlations [EOl Ell. 
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In addition to analytical techniques and experiments, 
several groups have utilized computer simulations of in- 
elastic hard spheres with both Molecular Dynamics 
H, [24| pBj a nd Direct Simulation Monte- 

Carlo ||28|, |29|, ^0|, 3l|, |32| techniques to investigate the 
statistical properties of granular gases. Using DSMC sim- 
ulations, Baldassarri et al. |^ have found velocity and 
density distributions that are qualitatively similar to our 
previous experimental results pT| . Recent work by Brey 
and Ruiz-Montero investigate how the second and 
fourth moments of the VDFs scale as a function of dis- 
tance from the driving wall, which until now, have not 
been experimentally tested. 

In this paper, we report on the statistical properties of 
a gas of inelastic particles constrained to two-dimensions. 
An inclined geometry reduces the gravitational accelera- 
tion acting on each particle which results in lower mean 
velocities. The combination of slow dynamics and high 
speed imaging allows us to accurately identify the parti- 
cle trajectories and collision events. By using velocities 
before and after a collision event, we measure the normal 
coefficient of restitution. We find that these quantities 
are found to be broadly distributed for the same impact 
parameters. By calculating the distance and time be- 
tween collision events we measure the distributions of 
free paths and times. We find that these distributions 
do not follow the result found from kinetic theory. The 
path and time distributions have an overpopulation of 
short distance and time bins, demonstrating the inher- 
ent clustering present in granular gases. We propose an 
empirical form that captures the distributions, which is 
then used to calculate the mean free path and free time 
as a function of density. The particle trajectories are also 
used to measure the mean square displacement, velocity 
auto-correlation, and diffusion rates. The distribution of 
particle velocities are measured with a variation in den- 
sity of an order of magnitude and show distinctly non- 
Gaussian behavior with no apparent universal form. We 
compare our results to recent experiments, as well as the- 
oretical and simulation treatments of equivalent systems. 

The paper has the following structure. In Section || we 
present the experimental apparatus and imaging meth- 
ods. Section [I] provides the overall system characteris- 
tics such as the density distributions and coefficients of 
restitution and inelasticities. We then present our analy- 
sis of the trajectories of the particles in Sec. IV. Finally, 
m Sec. we summarize our results in the context of 
granular kinetic theory and simulations. 



II. EXPERIMENTAL METHODS 

The experimental configuration [Fig. consists of a 
100 cr X 60 fj glass plane that is inclined at an angle /3 
with the horizontal. The particles are stainless steel with 
diameter a — 3.175 mm and a high degree of sphericity 
(6a/a = 10^*). The number of particles, measured in 
number of mono-layers Ni across the driving wall, is var- 
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FIG. 1: (a) Schematic diagram of the experimental setup. 
The inclined plane is a smooth glass surface, the side- walls and 
driving wall are stainless steel so that the particle-boundary 
collisions approximate those between particles. The driving is 
produced by a solenoid connected to the lowest side-wall. The 
angle of inclination /3, can be varied from j3 — — 8 degrees, 
the values of (3 we have chosen are 2° ± 0.1° and 4° ± 0.1°. 
(b) An image of the system taken from above. The bottom 
right corner is considered the origin of our coordinate system 
(0,0). The white bars allow us to track the position of the 
driving wall. 



ied between Ni — 1-5 in steps of one layer, (viz. from 
Np = 100 — 500 in steps of 100, where Np is the number 
of particles). The energy source is an oscillating side- 
wall, driven by a solenoid, that is located as shown in 
Fig. 0(a). The driving signal is a 10 Hz pulse with a 
velocity during each pulse of ~ 40 cm s^^. The driving 
frequency and amplitude were chosen to ensure that no 
phase dependence on the center of mass is observed (at 
frequencies below 2 Hz the particle positions are phase- 
locked with the driving). The signal is produced with a 
computer interfaced Aglient Technologies 33 120- A wave- 
form generator that is and subsequently amplified by an 
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FIG. 2: (a) Linked particle positions over 1365 time steps 
{Ni = 3). We can determine particle collision events with a 
high degree of accuracy from trajectories such as this, (b) 
The parabolic path of a particle. We use fixed values for Vx 
and g — 980 cm s~^ to measure g' = |(;sin/3. The fit gives 
f} = 2.2°, the deviation from the measured value of P is 9%. 



HP 6824A Amplifier. The inclination of the plane can be 
varied between (3 = 0° — 8° , for our experiments the angle 
was fixed at /3 = 2° ± 0.1° or 4° ± 0.1°. In the extreme 
case of /3 ^ 1°, the particles essentially cease to interact 
with the energy source and cluster at the side opposite 
of the driving. 

The particles are imaged using a Kodak MotionCorder 
SRIOOO high speed digital camera. We measure the posi- 
tions of all particles contained in the apparatus for 1365 
frames at 250 frames per second at full spatial resolu- 
tion of 512 X 480 pixels. These digital images are then 
transfered to a computer and analyzed using a centroid 
method that allows us to resolve each particle to sub- 
pixel accuracy. After each particle is located the positions 
are then connected in time to form continuous trajecto- 
ries for 5.46 s. Our coordinate system is such that the 
X, y axes are parallel and perpendicular to the driving 



respectively [see Fig. ||. A typical particle trajectory is 
shown in Fig. ||(a). Multiple collision events can be dis- 
tinguished with nearly straight paths between each event. 
A particle that freely rolls on the inclined plane will fol- 
low a parabolic trajectory [see Fig. |2|(b)]. The particle 
trajectory is given by 



y{x) 



5 X' 



(1) 



where, g is the acceleration due to gravity, Vx is measured 
from the width of the parabola, and the | factor is due 
to the moment of inertia for a solid sphere. 



III. SYSTEM CHARACTERISTICS 

A. Density Distributions 

The results presented will be given in terms of the num- 
ber of single layers across the cell, A^; and the angle of 
inclination /3, which the determine the area fraction (j) 
[see Table || . We measure (/) by defining a region of in- 
terest (ROI) that is centered about the peak in p{y), 
[Fig 11(b)] whose extent in the y-direction is limited to 
±10% of p{y). The ROI scheme excludes all particles 
that are within 3a of the side walls to ensure that clus- 
tering due to the side-walls does not affect our results 
[see Fig. ||(a)]. The over-plotted box in Fig. ^b) demon- 
strates the ROI definition for (j> — 0.13. A more stringent 
division of the system in the y-direction will be used when 
the behavior of the temperature, pressure and kurtosis 
are discussed in Sec The form of the density in the y- 
direction is similar to that found in Refs. ||lj, |3^ however, 
we find that the form of the tails of p{y) at higher values 
of Ni (3 deviate from Boltzmann distribution. This im- 
plies that the law of isothermal atmospheres breaks down 
for granular systems as we shall also see when we dis cuss 
the scaling of the granular temperature in Sec. IV E . 



Ni 



2.0 
2.0 
2.0 
2.0 
4.0 
4.0 



0.022 
0.068 
0.138 
0.191 
0.302 
0.581 



100 
200 
300 
400 
400 
500 



TABLE I: Experimental values of the number of layers A^i, 
the angle of inclination P, and the resulting measured value 
of (j). Np the number of particles in the system is given for 
clarity. 
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FIG. 3: (a) The density p{x) versus x for all iV;. The ob- 
vious clustering due to inelastic collisions at the side-walls is 
demonstrated here. Also, as A^; is increased the system be- 
comes more inhomogeneous in across the cell. This effect is 
most likely due to the onset of clustering instabilities that 
have been recently discussed H, 1^, ||. (b) The aerial den- 
sity plots p{y) for each Ni and /3 on a log- linear plot. The 
density </> is measured in a particular area by integrating p{y) 
over that region of interest. The total are under each curve 
corresponds to the average area fraction for that particular 
Ni . The solid line shows an exponential fit over the tail of the 
distribution of Ni = 200. However, we will demonstrate that 
the isothermal atmosphere is not obeyed for any density. 



B. Particle Collisions 

We identify collision events from the trajectories by 
using the following algorithm. Velocities are constructed 
as finite differences Vj = where Ax = x(tj) — x(ii) 
and the subscripts i,j represent positions separated by 
the time difference At — 4 ms. All velocity vectors are 
compared sequentially to find direction changes given by 



where v = v/|v| the unit vector of the calculated veloc- 
ity. If 20° < t/j ^ 180° the proximity of all particles at 
the same time instant is checked. If a particle is found 
within a radius a + Act, whose velocity also satisfy Eq. 
(^) it is considered as a candidate for a collision. To as- 
sure that re-collisions are not occurring, we maintain a 
record of the identity of the previous collision partner. 
We then ensure that those particles can re-collide if and 
only if the partner particle has undergone a collision with 
yet a third particle. If particles pass these requirements 
then a collision has occurred. To extend the algorithm to 
include collisions with the boundary walls we first check 
if Eq. (^) is satisfied. We then check if the particle's 
center is within a + Aa of a boundary and it's velocity 
component perpendicular to the wall is reversed. 



C. Coefficient of Restitution 

The loss of energy in a collision is given by the co- 
efficient of restitution. Two particles that undergo an 
inelastic collision with a relative velocity between parti- 
cles Vi2 = vi — V2 , will obey the reflection law v^j • = 
—a Vi2 • (T, where a is the normal component of the resti- 
tution coefficient and a is the unit vector connecting the 
centers of the particles. 

Having an efficient method for collision identification, 
we are able to measure the relative velocities of two par- 
ticles before and after collision events. The coefficient of 
normal restitution during a binary collision is given by 



(Vl2 • 0-) ' 



(3) 



where the over-bar denotes average over three pre/post- 
collisional velocities measured in the ROI described above 
[see Fig. ||(b)]. The angle between the relative velocities 
of two colliding particles is given by 



= cos 



12 



'12 



(4) 



Thus we can characterize the coefficient of restitution as 
a function of 9. The probability distributions P(q;) for 
60° < 6* < 180° for each Ni (3, are shown in Fig. |(a- 
f). Data for < 60° suffers from a lack of statistics 
and therefore is not included. Each graph represents the 
probability of the inelasticity having a value a for a range 
of 6* + A0, where A9 = 2°. P{a) follows a very broad 
distribution of values over all 9, and have a decreasing 
mean value as function of 6 [see Fig. |(a)]. Thus we 
find that the coefficient of restitution can have a broad 
distribution of values for the same impact angle. 

We also measured the energy loss due to a collision as 
a function of Ni (3. The ratio of the magnitudes of the 
relative velocities before and after a collision. 
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FIG. 4: The distribution of the normal component of resti- 
tution a versus 60° < 9 < 180°, the relative angle of inci- 
dence between particle velocities, (a) A*'; = 1 ,P — 2.0, (b) 
Ni = 2,p = 2.0, (c) iV; = 3 , /3 = 2.0, (d) Ni = 4 , /3 ^ 2.0, (e) 
Ni =4,(3 = 4.0, (i) Ni = 5(3 = 4.0. The value of the z-axis 
for each graph is the probability of a collision giving a value 
of a in a range of 6 + A6, where A9 = 2° 
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FIG. 5: (a) The mean value of the distributions of a shown 
in Fig. ^ averaged over 60° > 9 > 180°, as a function of the 
average covering fraction <j). The bars indicate the spread in 
the distribution, (b) The distribution of energy inelasticities 
given by Eq. (|) for {>)Ni = 1,(3 = 2.0, {<i)Ni = 2,(3 = 2.0, 
{A)Ni =3,(3 = 2.0, {o)Ni =4,(3 = 2.0, (□)iVi = 4,p = 4.0, 
(o)A^; = 5 (} = 4.0. Each distribution is shifted vertically for 
clarity. 



determines the energy restitution coefficient, (ry^ = if 
all 9 are averaged). Figure ||(b) shows the distributions 
of measured values of rj shifted for clarity. We find that 
a peak exists at a value that is consistent with a^. Fur- 
thermore there exists a power law tail for values of 77 > 1 , 
which has been interpreted as a random inelasticity ^T\ . 
The appearance of a tail at high rj implies that the rota- 
tional degrees of freedom are actively transferring energy 
to translational motion during a collision. 

IV. RESULTS 

A. Distributions of Paths and Times 

We measure the distribution of paths lengths from the 
the geometric distance between collision events defined 
in our ROI at each Ni P [Fig. ^a-f)]. By basic kinetic 
theory arguments pq|, the distribution of path lengths 



for an elastic hard-sphere gas (and by a similar treatment 
the distribution of free times) is given by, 

P(/) = (2%/2 0) e-^v^"^'. (6) 

The distribution therefore should follow a simple expo- 
nential form depending only on the density. However it 
is clear from the dashed lines in Fig. ||(a-f) that the sim- 
ple form given by Eq. (|6|) does not describe the behavior 
over all I. 

The distributions of times between collisions P(t) 
[Fig. ^(a-f)] is also measured and shows similar behav- 
ior to that of the path length distributions, that is an 
overpopulation of the short time bins. This should be 
expected from the simple relationship between the dis- 
placement and the time. However, it is noteworthy to 
mention that the ratio of Z/r versus the path length I, 
is not a constant over all values of I, implying that the 
average speed of the system depends on the distance or 
time between collisions. 
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FIG. 6: The probability distributions of path lengths P{1) 
versus on a log- linear scale, and inset log- log scale (a) Ni — 
1,13 = 2.0, (b) Ni = 2,13 = 2.0, (c) TV, = 3 , f3 = 2.0, (d) 
Ni = 4,/3 = 2.0, (e) N, = 4:,f3 = 4.0, (i) N, = 5 , f3 = 4.0. 
The dashed line shows the theoretical form given by Eq. (^) 
derived for elastic particles, and the solid line is an empirical 
fit given by Eq. (^). Table |l| shows the fit parameters. 



Elastic hard spheres will have a mean free path that is 
simply I = vf, where v and f are the average speed and 
collision time respectively. Also, the mean free path can 
be derived directly from the distribution of path lengths, 
I = I P{1) dl, where P(Z) is given by Eq. g[ Grossman 
et al. ^ have interpolated how the mean free path for a 
granular system should be modified to account for higher 
collision rates due to increased density. Although the in- 
terpolation gives a qualitatively accurate correction for 
passing between the high and low density limits, the ac- 
tual distribution of path lengths has not been measured 
or calculated for a granular gas. 

We have found an empirical form that well describes 
the measured distributions of path lengths and free times. 



Pil) = a{l) 
P{t) = a{T)-^e 



(7a) 
(7b) 



where a, b^c for the path lengths and free times are shown 
in Table y for all Ni (3. This form appears to capture 
both the short / and r power-law behavior. In the dilute 
case the form returns to the theoretical prediction for 
larger path lengths. 



FIG. 7: The probability distributions of free times P(r) versus 
r, on a log-linear scale {&) Ni = 1 , 13 = 2.0, (b) iV; = 2 , /3 = 
2.0, (c) iV, =3,/3 = 2.0, (d) iV, =4,/? = 2.0, (e) Ni^A,l3 = 
4.0, {I) Ni ^'^,(3 = 4.0. The solid line is a fit given by Eq. 
(ib) 



From the distribution of path lengths and free times, 
we calculate the mean free path and time by utilizing the 
fitting form and it's parameters. The ratio of the mean 
free path to the mean collision time should determine the 
average speed v in the ROI where the distributions are 
measured. We have taken the ratios of the integrated 
distributions, 



I _ j^lP{l)dl 



(8) 



f j^TP{r)dr' 

and compared that to the average of the speed distribu- 
tion {v^^y) in the same ROI. Figure [s] the both measure- 
ments for all Ni (3. The agreement is within 10% over the 
entire range of iV, (3 indicating that the proposed forms 
in Eqs. (0a, b) quantitatively capture the behavior of the 
distributions. 



B. Velocity Autocorrelation 

The velocity autocorrelation function (VAF) is com- 
puted for the X components of the velocities within an 
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FIG. 8: The average speed measured for each A^; fi. (o) v 
obtained from Eq. (p|), and (o) {v) measured from the mean 
of the speed distribution. The two independent measures give 
similar values over the entire density range. 



ROI by using the following [B9|: 



N N i 



jj=0 At=l 



E 



,(to) •v,j(to + At) (9) 



where, Ng is the total number of data sets, Np is the 
number of particles and tmax is the total number of time 
origins. Figure ^(a), shows the measured values of the 
VAF normalized by (v(0)^) in our system. 

In simulations of hard sphere fluids, Alder and Wain- 
wright 1^ first found that the form for the VAF was 
strongly depended on the density of the system. For very 
low densities the the characteristic form of the correlation 
function was given simply by 



, -t/TE 



(10) 



where te is the Enskog collision time. If the density of 
the system is increased however, the form of Eq. ( p^ ) 
breaks down and Cy{t) can become negative with long 
range tails due to the caging of particles by their neigh- 
bors. We find, that the lowest density case becomes, 
and remains negatively correlated after the decay from 
(v(0)^) [Fig p|(a)]. This appears to be in contradiction 
with Ref. ||4C|] but is due to a finite size effect. That is, 
the particles are interacting frequently with the side- walls 
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FIG. 9: (a) The velocity auto correlation function Cvit). (b) 
The mean square displacement in the x-direction for each A*'; (3 
for t = Q — 2.5 s. Inset: The short time behavior indicated 
by the box in the main figure, (c) The diffusion constants 
calculated for each cf>, where (o) corresponds to the numerical 
integration of Cv from (a), and (o) corresponds to the least 
squares fit of C^2 from (b) . The dashed line shows the kinetic 
theory result for a fixed temperature, which is given by the 
average measured value over this range of (j>. 



Ni 




a{a) 


m 


c(c) 


1 


2 


0.031(0.0027) 


0.428(0.511) 


0.154(2.969) 


2 


2 


0.025(0.0025) 


0.603(0.665) 


0.393(6.024) 


3 


2 


0.008(0.0003) 


0.932(1.203) 


0.742(8.306) 


4 


2 


0.003(0.0008) 


1.258(1.209) 


1.489(16.91) 


4 


4 


0.002(0.0002) 


0.970(1.043) 


3.171(26.45) 


5 


4 


0.002(0.0005) 


1.096(1.384) 


2.887(28.33) 



TABLE II: Fitting parameters for Eqs. (ya,b). The values are 
arranged a{a) for P(Z)(P(r)), respectively. The (■ ■ ■ ) corre- 
spond to the values for P{t) 



at low densities, which reverse the sign of velocity vec- 
tors, thus leading to the observed anti-correlation. The 
predominance of the sidewall interactions are screened for 
the intermediate densities due to the increased number of 
particle-particle collisions, therefore no anti-correlations 
are observed. 



C. Mean Square Displacement 

To determine the mean square displacement of the x- 
component of the particle positions [Fig ^(b)] , by using 
the sequential time zeros for each trajectory [p9[ given by 
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E 

=0 At=l 



E 



|Xy(io)-Xy(to + Ai)P. 

(11) 

Where Ng is the total number of data sets, Np is the 
number of particles and tmax is the total number of time 
origins. For these measurements, the ROI is allowed to 
increases in size along the y-direction as Ni decreases. We 
have chosen to make this increase to ensure that particles 
at low Ni have had the opportunity to undergo a collision 
while under consideration. 

The long time behavior of Cx^ for each Ni (3 [Fig. ||(b)], 
displays linear dependence on time, indicating diffusive 
behavior. However for iV; = 1, Cx^ clearly shows a cross- 
ing from one linear regime to another, which may be a 
possible indication of finite system size for low density. 
For short times, [Fig. ^(b)/nsei] the behavior is ballistic 
as indicated by the quadratic increase of C^2 in time. 
As Ni is increased, the range of the ballistic regime dra- 
matically decreases indicating a decrease in the Enskog 
collision time te- The ballistic and diffusive regimes are 
consistent with what is expected for kinetic theory of 
elastic, finite-sized particles. 



D. Self-Diffusion 

The self diffusion constant D, can be determined for a 
system of particles by either evaluating the time integral 
of the velocity autocorrelation function, 



D 



(12) 



or using the relationship between the mean square dis- 
placement of the particles and time over long times. 



(13) 



where d is spatial dimension. From kinetic theory [ p8[ , 
the diffusion constant of a two dimensional gas is calcu- 
lated as: 



D = 



ttT\ 

805(a) y^J 



1/2 



(14) 



where, g{a) is the radial correlation function at con- 
tact [M given by. 



lir = cr) 



16-70 
16(1 - 0)2 



(15) 



By numerically integrating the curves in Fig. ||(a), and 
performing a least squares fit to the data in Fig. ||(b) af- 
ter the ballistic regime, we obtain the self-diffusion con- 
stant [see Fig. ^(c)]. We find that the values for the 
self-diffusion from Eqs. (p^JT^) are self consistent. The 



solid line in Fig. g(c) shows the form of Eq. (|lj) with the 
temperature T given by the granular temperature. The 
granular temperature is defined by 



(16) 



denote 



where, m is the mass of the particles and ( 
averages over the component distributions [see Sec. [VE|. 
The other constants in Eq. |l^ are determined from sys- 
tem parameters. The theory for the diffusion of elastic 
particles, given by Eq. ([l^ ) closely matches our results 
all (p. Thus we show that the effects of inelasticity on the 
self-diffusion are small. 



E. Velocity Distributions 

The distribution of the x- and y-components of the 
particle velocities are plotted in Figs. [l0|-[l2|(a-f). The 
distributions correspond to velocities that are measured 
within a region of interest. The ROI is defined by making 
a narrow slice across the y-direction that is centered upon 
the peak in p(y) while excluding particles lying within a 
distance of 3<j from the side walls. We utilize this ROI 
to ensure that large gradients in p{y), and the clustering 
produced by the side-walls, do not affect the measured 
VDFs. Each distribution correspond to 2 x 10^ unique 
velocities that are found within our ROI. 

The velocities of elastic particles follow a distribution 
given by the Maxwell-Boltzmann form. 



P(v) = {27:kBT) 



-<i/2g-vV2fc£ 



(17) 



where, d is the dimensionality of the system and T is 
the temperature of the heat bath that the system is in 
contact with. Hence, if a system of particles is at equi- 
librium, its temperature determined by the width of the 
distribution of particle velocities. Equation ( |l7| ) is fit to 
the data for the a;-components of the velocities, and is 
shown on both linear and logarithmic scales [Fig. |lO| , pl| . 
We observe that the form given by Eq. ( |l7|) displays de- 
viations both at low and high velocities. The distribu- 
tions of velocities are normally displayed in a log-linear 
fashion to accentuate the tails of the VDF, however this 
suppresses the deviations at low velocities. By plotting 
the distributions on a linear scale we display the more 
statistically significant deviations from Eq. (p7|). 

In a recent experimental work , a two dimensional 
collection of particles is driven into a steady state. Using 
analysis techniques that are similar to ours, the authors 
proposed a governing form for the VDF given by 



i?(v)=Ae-^l^-/^^l" 



(18) 



where A and B are constants and Tx is the a;-component 
of the granular temperature defined in Eq. (p6|). They 
claim to have seen a universal VDF which they parame- 
terized by a single value, regardless of the system density 
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FIG. 10: The velocity distribution functions P{vx) versus Vx 
on a linear-linear scale, (a.) Ni = 1 , f3 — 2.0, (b) Ni = 2 , f3 — 
2.0, (c) iV, = 3,13 = 2.0, (d) Ni = i,(3 = 2.0, (e) iV, = 
4 , /3 = 4.0, (f ) iV; = 5 , /3 = 4.0. The solid curves are a least 
squares fit to a Gaussian form given by Eq. (0). Note that 
the deviation from a Gaussian distribution extends all the 
way to the lowest velocity bins. Each distribution correspond 
to ~ 2 X lO'' unique velocities that are found within the ROI 
defined in the text. 



FIG. 11: The velocity distribution functions P{vx) versus Vx 
on a log-linear scale.(a) iV, = 1 , /? = 2.0, (b) iV; = 2 , /3 = 2.0, 
(c) Ni =3,13 = 2.0, (d) Ni=4,f3 = 2.0, (e) iV, =4,/? = 4.0, 
(i) Ni = 5 , f3 = 4.0. The solid curves are a least squares 
fit to a Gaussian form given by Eq. (|l^). Here the apparent 
deviation in the tails of the distribution functions are present. 



or the value of the inelasticity of the particles. From our 
VDFs, whose corresponding densities range over an order 
of magnitude and where the average inelasticity varies by 
nearly a factor of two, we cannot find any single param- 
eter fit that describes the overall form. 

The VDFs for the y components, -P(vj,) versus y^, for 
each TV; P in our ROI are also measured [see Fig. |l2(a- 
f)]. The VDFs are highly skewed by the asymmetry in 
the driving against the direction of gravity. To identify 
the effects that the asymmetry in P{'Vy) has upon P(vj;), 
we have separated the v-c distributions by the sign of Vj^, 
i.e. P{yx\ + Vj,;— Vj,). We have found that the form 
for these conditional distributions are not affected by the 
sign of Vy, however we do note that their widths differ. 



with 



We also measure the a;-component of the granular tem- 
perature [see Eq. (p^], to probe the scaling behavior 
of the velocity distributions. Figure |l^(a) shows the mea- 
sured granular temperature as a function of distance from 
the driving wall. At low densities {Ni < 2), Tx{y) initially 
increases and then decays. In contrast, for {Ni > 3), 



Tx{y) has a distinct minimum. We note that Tx{y) never 
reaches a constant value and the minimum (maximum) 
does not correspond to the peak in p{y). 

To further show the non-universality of the VDFs, we 
plot the kurtosis as a function of distance from the driv- 
ing wall. The kurtosis is obtained by the following: 



(19) 



ivl?' 



If the velocity distribution is a Gaussian then 7 = 3, 
shown by the solid line in Figure 13(b), and if the dis- 
tribution is given by Eq. ( [l8| ) then 7 = 3.576. We find 
that the the measured values for 7 exceed the value for 
a Gaussian and also vary as a function of distance from 
the driving. This analysis is consistent with our previ- 
ous results [^ and recent MD simulations of Brey and 
Ruiz-Montero p7| that closely mimic our experiment. 
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FIG. 12: The velocity distribution functions P(vy) versus Vy 
on a log-linear scale. (a) iV, = 1 , /3 = 2.0, (b) A^j = 2 , /3 = 2.0, 
(c) Af, =3,/3 = 2.0, (d) TV, =4, /3 = 2.0, (e) iV, = 4 , /? = 4.0, 
(f ) A^! = 5 , /3 = 4.0. The large skewness in the distributions 
for the negative values of Vy is due primarily to the driving 
from the bottom wall. Particles that are moving in the —y 
direction are leaving the moving wall. 



F. Equation of State 

The equation of state for ideal gases relates the pres- 
sure to the temperature and the density, 



(20) 



where n is the number density and fc^ is Boltzmann's 
constant. If we assume that kinetic theory is valid for 
a granular gas, we can immediately relate the average 
squared speed, (u^) of the particles to the temperature. 



(21) 



for each degree of freedom. Figure ^(a) shows the pres- 
sure, P = ^{Vy) as a function of distance from the driv- 
ing wall. Due to the effects of gravity on the particles, the 
density should follow the well known atmospheric law. 



(22) 



which assumes a constant temperature. We find that the 
temperature is not constant for any Ni (3 [see Fig. ^ . 
This is also consistent with our observations of the den- 



sity distributions in Sec. [II A where p{y) deviates from 
the form of Eq. (151 




FIG. 13: (a) The granular temperature \m{y%) as a func- 
tion of distance from the driving wall for each Ni f3. If the 
isothermal atmosphere condition was satisfied these would be 
constant values for all y above the peak in p{y). For values of 
Ni > 2 the temperatures follow a non-monotonic form that 
has a distinct minimum, (b) The kurtosis, 7^ measured from 
P{'Vx) as a function of the distance from the driving. The val- 
ues give n by a Gaussian (solid line) and the form proposed in 
Eq. (|l8[) (dashed) are only attained very far from the energy 



Momentum balance implies that the gradient of the 
pressure is related to n by the following equation: 



dP 
dy 



~nmg 



(23) 



where, m is the mass of a particle and g is the acceler- 
ation of gravity. Due to the non-vanishing gradient of 
temperature the general form of the pressure gradient 
must be taken into account 



dP dn 
dy "dy 



dTy 
dy 



(24) 



We find that Eq. (|2J) is indistinguishable from the nu- 
merical derivative of d{nTy)/dy. 

We have measured the pressure gradient acting on a 
particle held at a particular y by evaluating. 



Ty \ d jnTy) 
dy 



y_ 

nm 
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(25) 



where g' — |(7 sin(/3). Figure ^(b) shows the LHS of Eq. 
(H) for /3 = 2, and Fig. |l|(c) for /3 = 4. The solid lines 
correspond to the values of g' with /3 = 2, 4 found by 
utihzing Eq. (|I]) and the data in Fig. ||(b). We find that 
the measured values systematically overestimate the ac- 
tual values for g' in the region where the density reaches 
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FIG. 14: (a) The pressure P{y) = pT as a function of y, the 
distance from the driving wall, (b) The ratio of the mass 
density to the granular pressure force — 4^ as a function of 

^ ^ mp ay 

distance from the driving wall. The solid line corresponds to 
g' for /3 = 2° (c) The same as (b) for /3 = 4°. The obvious 
deviation for regions of high density show a breakdown of the 
simple treatment of the granular equation of state in regions 
of high density. 

it's maximum. Our interpretation assumes a dilute gas, 
therefore the deviations near the peak in p{y) are not 
surprising. 

To incorporate the effects due to increased density, we 
have obtained the pressure from the interpolated equa- 
tion of state derived by Grossman et al. |6[ : 



Uc — n 

where, Uc is the close packing number density. Utilizing 
Eq. to calculate the gradient of the pressure [Eq. (p3|)] . 
However, we find that the disagreement persists between 
the LHS of Eq. ( [25[ ) and the measured value of g' near 
the peak in p{y). 

V. SUMMARY AND CONCLUSION 

In this paper we have presented a statistical analysis 
of an inelastic gas that is constrained to two dimensions. 



Utilizing high speed digital image processing we perform 
long time tracking over a broad range of densities. Not 
surprisingly, we observe that the statistical properties of 
inelastic gases deviate from expectations of the kinetic 
theory for smooth elastic particles. The most apparent 
discrepancies arc found in the distribution of free paths 
and times and the distribution of particle velocities. 

To characterize our system we measure the effective 
coefficient of restitution from the relative pre- and post- 
collision velocities of particles undergoing binary colli- 
sions. We find that the normal component of restitution 
and the energy inelasticity are not single- valued, but have 
a distribution of values even for the same impact param- 
eters. The mean value of the normal components of resti- 
tution systematically decreases with the system density. 
We also find that the energy inelasticity can take on val- 
ues greater than unity, demonstrating a transference of 
energy from the rotational to linear degrees of freedom. 
In a recent numerical work Barrat and Trizac have mea- 
sured the projected one-dimensional coefficient of resti- 
tution Ref. Il^, Q and the energy inelasticity. Their in- 
terpretation is that the coefficient of restitution and the 
energy restitution are random variables that characterize 
collisions, consistent with our findings. 

The distribution of path lengths and free times are 
shown to have an overpopulation of the short distance 
and time bins. We have proposed an empirical form in 
Equations (j^ a,b) which capture the overall behavior of 
the observed distributions of paths and times. By inte- 
grating these distribution functions, we are able to mea- 
sure the mean free path and mean time. The average 
speed obtained from the speed distribution and from the 
mean free path and time arc in close agreement. Inspired 
by these finding, Paolotti et al. ||2^ have reported similar 
results for the mean free time in a simulation that mimics 
our system. 

Particle diffusion constants are measured from two in- 
dependent long time averaged correlation functions. The 
mean square displacement and the velocity autocorrela- 
tion function are calculated. By then performing least 
squares fitting and numerical integration to these quan- 
tities respectively, the self diffusion over a broad range 
in density is calculated. We find that the diffusion con- 
stants are similar to that of a two dimensional gas over 
this density regime. Therefore, long-time averaged corre- 
lation functions seem to accurately capture the diffusive 
properties of granular gases. 

We find that the distribution of particle velocities per- 
pendicular to the direction of driving does not have a 
universal form, but depends on both the density and the 
inelasticity. In addition, we find a distinct asymmetry in 
the VDFs parallel to the driving direction. We measure 
the granular temperature as a function of distance from 
the driving source and find non-monotonic behavior. For 
low densities, the granular temperature has a distinct 
maximum and for high densities there exists a distinct 
minimum. The temperature inversion at higher densities 
has recently been described via granular hydrodynamics 
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by Ramirez and Soto [Q. However, the crossover from 
T{y) having a maximum (for low densities), to T{y) with 
a minimum at high densities has not been discussed in 
any kinetic or hydrodynamic models. 

By using kinetic theory and simple hydrodynamics we 
have tested the force balance between the gradient of the 
pressure exerted by a granular gas on a particle and the 
force due to gravity [Eq. (|2^)] . We find strong deviations 



in the regions of high density. A simple hydrodynamic 
form, that describes the behavior over all densities, is not 
yet available. 
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